Double Scaling Limits and Catastrophes of the zerodimensional 
O(A^) Vector Sigma Model: The A-Series 



W. Riihl* 

Fachbereich Physik, Universitat Kaiserslautern, 
D-67653 Kaiserslautern, Germany 

J^TgX- version of preprint KL-TH-94/21 
September 1994 



Abstract 

We evaluate the partition functions in the neighbourhood of catastrophes by saddle point 
integration and express them in terms of generalized Airy functions. 

1 Introduction 

The vector models have attracted interest as one-dimensional quantum gravity theories just as 
the matrix models are interpreted as representing two dimensional quantum gravity theories. A 
connection with polymer models has also been emphasized from the outset Their double 
scaling limit has been studied with the usual renormalization group and 1/N expansion methods. 
In the zero dimensional case the beta function and the free energy have been calculated for large 
exactly. 

We will show that the double scaling limits for this most elementary model can be calculated 

1 

exactly, i.e. asymptotically to any order in a A "i+i -expansion. Catastrophes are singularities 
of differentiable maps |^ and by diffeomorphisms can be transformed to canonical forms. We 
will study such canonical forms only. It does not make sense to reshape these canonical forms 
by application of diffeomorphisms. 

There are elementary and nonelementary catastrophes. The elementary ones are ordered 
into A, D and E cases Whereas vector models with one vector field can exhibit only A- 
series catastrophes {Am, rn E iV), models with two vector fields can also posses D or E series 
catastrophes . The nonelementary catastrophes show up also first in two-field models. 

Application of catastrophe theory to zero dimensional 0{N) sigma models leads to a wealth 
of useful information which can be used as a guideline for studies of more complicated models. 
The diffeomorphisms which are basic in catastrophe theory replace the 'reparametrisations in 
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coupling constant space' ; those properties of catastrophes which hold true independently of 
diffeomorphisms (characterize cosets of the diffeomorphism group) map onto universal features 
of phase transitions. This is of course not new. Nevertheless our analysis will lead to some new 
insights. 

We define the zero dimensional O(A^) vector models by 

2N{g) = J exp I - iV5(<I))|, : Lesbeque measure, ^ G ]R]n (1.1) 

where g{^) is O(A^) invariant and has the asymptotic expansion 

1 oo , 

Constraints on the {gk} lead to catastrophes that dominate the large behaviour of the partition 
function (|1.1| ) through saddle point expansions. Main ingredients in these expansions are (real 
nonoscillating) generalized Airy functions. 
The free energy is defined by 

FN{g) = -j^logZ^ig) - \log^ (1.3) 



where this normalization is such that 



FN{g) 



= (1.4) 

gk=0 Vfc>2 



The simplest catastrophe is of Ai type (or Morse or Gaussian). We set 

51 = 1, 52 > 0, 5fc = 0V/c>3 (1.5) 
then Zfyf can be expanded into a 1/A^ expansion around a Gaussian saddle point 



1. 1 + vA 1 1 1 , , 

F„fe)^_ = ^log— + (1.0) 

with 

A = 1 + 4^2 (1.7) 
This function satisfies the (large A^) Callan-Symanzik equation. 

{^-^-f^(s2)^^+i{g2)^FN{g2) = RN{g2) (1.8) 

Assume we know that 

7(52) = 1 (1.9) 
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for the free energy. We can continue F]\f{g2] 
neighbourhood of 



N=oo 



analytically off the positive real axis till the 



1 

1 



(1.10) 



where F]\f{g2 



N=oo 



has a branch point in the variable A 



FN{g2 



N=oo 

HO) 



g{A) + A-2 hiA) 
1 

~3 



(1.11) 



both functions g, h being analytic. 

Now the singular part is defined to satisfy the homogeneous Callan-Symanzik (renormaliza- 
tion group or RG) equation. It follows 

d 



and 



(3(92)-' 
(3(92) 



Rn{92 



log A2/i(A) 



■^{92-9*) + 0({g2- g.f 



(-te)|^ + l).(A) 



The results on (3, Fn, Rn given in Q are thus reproduced. 

In the neighbourhood g2 — g* it is guessed that the singular part of -F/v(52 
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= cxD is equal A2/i(A)) possesses the series expansion 



'smg 



FN{g2 



'smg 



J2 ahN~^'{g2-g.\ 



2-70-71^ 



h=0 



This is based on the fact that each term of ( |1.15 ) satisfies the RG equation when 



71 



/3'(S*) 



3 

2' 



70 



7(g«) 



(1.12) 
(1.13) 

(1.14) 
(which at 

(1.15) 
(1.16) 



However: this argument is too simple, since the RG equation does not determine h to be an 
integer. If we apply derivation of h to any term in ( 1.15D we obtain also a permitted contribution 



a'h N-^ (52 - (7*)2-^o-^i'^ log \N{g2 - g.V' 



(1.17) 



In fact such term for h = 1 exists. 

The g2 parameter defines a curve of Ai catastrophes which at g2 = g* ends in a A2 catastro- 
phe. But at this catastrophe and in a neighbourhood we can derive -F/v(5'2)sing directly (all of 
Fn{92) in fact). The result (Section 3) is 



-PV (52) sing 
C 



^logBi(C) 
{92 - 9*) 



(1.18) 
(1.19) 
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and Bi is the usual Airy function. Moreover the expansion (1.15) is asymptotic only and 



a'l = +1 (1.20) 

In section 4 we deform A2 catastrophes by 53. The aim is to render the partition function for 
canonical A2 convergent. A similiar problem with convergence appears for all Am, m even, and 
can be solved the same way. In section 5 we discuss the canonical Am catastrophes for general 
m and in section 6 we derive the corresponding /^-functions. 

2 A pedagogical example of mathematical interest 

The confluent hypergeometric function 

is a well studied transcendental function known in mathematical physics since more than hundred 
years. Nevertheless the asymptotic behaviour in the limit 

a — > 00, z ^ 00, 7 fixed (2-2) 
f = (2.3) 

1 2 
^ —-, so that (1 + 4^)z3 is again fixed (2.4) 

is not covered in the textbook literature (see e.g. Luke's otherwise extremly useful treatise 
This is a typical 'double scaling limit'. The standard approach is a saddle point integration 
technique. 

The integral representation for iFi is (|Q|, 9.211.2) 

iFi{a;r,z) = T^r^^ f' dtt-\l-ty-'exp{zt + alogt-alog{l-t)} (2.5) 

r(a)r(7 - a) Jo 1^ J 

for 5R(a) > 0, 3?(7 — a) > 0. Though the asymptotic region ( |2.2D lies outside the convergence 
domain of (|2.5| ), the saddle point expansion we shall derive remains valid. We consider only 
the contribution of the saddle point and not those of the boundaries. In physical applications 
one must be more careful. The relevant contribution is always the one which dominates the 
asymptotic behaviour. 
Define using ( p. 3D 

fit) = f + eiog^ (2.6) 

as 'phase function'. There are two extrema of f{t) at t± if ^ > — 1/4 and none if ^ < —1/4 

t± = (2.7) 
A = 1 + 4^ (2.8) 
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For A = we obtain a point of inflexion at 

io = ^ (2.9) 

We expand f{t) around Iq 

fit) = f{to) + {t-to)f'{to) + ^{t-toff"'{to) + 0[{t-to)') 

with 

f{to) = \. f'ito) = A, /"'(to) = 32e 
where we may approximate 

nto) = -8 

Now we scale the integration variable so that 

t — to = 

A = {Azyi 

Thus the leading part of iFi coming from this saddle point is 
where <1*(C) is of greatest interest for us 



Jc 



_ 1 2 

C = 4 3Z3A 

As a rule we use the average over two 'least deformed real axis' contours on which 
converges. Define 

r £ Q : Cr = contour from zero to e^'^*''oo along a ray 



Then in (2.17) we set 



With the ray integrals 



C = Co-l{Ci+C_i] 

2 3 3 



JCr 
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we get 

m = R^^\c)-UR?\0 + R'^'\iC)) (2.22) 

Z ^ 3 3 ^ 

Using small C expansions to identify functions we find ([||, 10.4.3) 

$(0 = vrBi(C) (2.23) 
(C)-i?i'UC)) = vrAi(C) (2.24) 

From ([^, figs 10.6 and 10.7) we see that 

Bi(Co) =0, Co = -1.173 (2.25) 

and Bi((") oscillates for C < Co and is positive for C > Co- The function Ai(C) oscillates everywhere. 
This justifies the choice of of the contour C ( |2.20| ) for C > Co. 

The Airy functions possesses a large C asymptotic expansion (|Q, 10.4.18 and 10.4.63) which 
is obtained by keeping —\ fixed and by expanding /(t) around the extrema t± (|2.7D . This is 
a saddle point expansion of A\ type and can be used to obtain information on the domain where 
$ is real non oscillating. Moreover we need it to recover the expansion ( 1.15 ) and to correct it 
(Section 3). 

The residue i?(t) of /(t) in ( ^.10| ) which is 0{(t — to)^) ai^d has been neglected still can be 
expanded systematically 

oo oo 

exp{zi?(0} ^ 1 + ^ ^ a„,fcz"(t-to)''^ (2.26) 

n=l k=2n 

with ttn^k polynomials in the derivatives at to. Moreover the function 

B{t) = t-^{l-ty-^ (2.27) 

can be expanded 

oo 

B{t) = Bito){l + Y.brit-toY) (2.28) 

Inserting both expansions into the partition function and submitting it to the same procedure 
as before we get ( ^.16 ) with $(C) replaced by 

oo oo oo 

E E E 2"§(2^+^)a„,fc6,z"-^(2fc+.)^(2fc+r)(^)^ (^^^^ ^ ^ 0) (2.29) 

n=0 r=0 k=2n 

which is an asymptotic expansion in powers of z~3 . 
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3 Elementary catastrophes , in particular A2 



Families of polynomials 

Cit + C2i2 + . . . + Cfc-it'-' ± (3 1) 

define a catastrophe of type with {(i, (^2; • • • > Cfc-i} as deformation parameters. Saddlepoint 
expansions around deformed catastrophes are dealt with in the encyclopadic treatise Q. Many 
questions relevant to our problem remain unanswered. In particular we would like to know 
where the generalized Airy functions $(Ci; C2, ■ ■ ■ , Cfc-i) are real positive. This domain in M^-i 
is bounded by a (A; — 2)-dimensional surface which can only be determined numerically. In 
section 5 we will learn that this question is relevant for even k only. 

In this connection the asymptotic behaviour of the functions <I> is of interest. But even for 
the Pearcy function Q 

^(Cl,C2) 

which is the last one in this series carrying a name, the above series are unknown. 
Let us return to ( |1.1[) now and perform the angle integration 

N 

with 



f{t) = logt-J2fi' (3-3) 

k=l 



We shall always normalize 



gi = 1 (3.4) 

We want to evaluate ( p.2p in the limit of large A^. If 

f'ito) = 0, /"(to) / (3.5) 

we have a Gaussian integral as leading term (^1 catastrophe or Morse singularity) implying a 
pure expansion. The A2 case arises if 

f'ito) = 0, /"(to) = 0, /"'(to) ^ (3.6) 

For this case to occur it is sufficient to have g2 as only coupling constant 

gk = 0, k>3 (3.7) 



The integral ( |3.2D can then be evaluated as a sum of two iFi-functions, which can be treated as 
in the preceding section. But the result thus obtained can be directly derived from the integral 
(3.2) by a saddle point technique. 
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Solving Q with Q gives 



to = (-52) ^ = ^ (52 < 0) 
/"'(to) = 16 

The deformation of the catastrophe is achieved by one free parameter, say A ( |1.7| ) 

A = 1 + 452 

/' = (1-A)^ - 1 

whereas 

to = (-52)"^ 
still holds. If we expand f{t) around to ^-t A ^ 

1 8 

fit) = f{to)- -A{t-to) + -{t-tof + lemamdei 



and scale the integration variable 



t — to = Xt] 
X = to 



we obtain as leading part of the partition function 



2sk 



N 



JV 



2^k 



where 



and 



$(C) = J drie~''''+^ 
c 

The contour C must be chosen such that under replacement of 
becomes identical with ( |2.24| ). 



(3.8) 
(3.9) 

(3.10) 
(3.11) 

(3.12) 

(3.13) 

(3.14) 
(3.15) 

(3.16) 

(3.17) 
(3.18) 
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In the case of the Airy functions asymptotic expansions for large C need not be obtained 
from a saddle point expansion but Q, 10.4.63 tells us that 

oo 

Bi(C) ^ n-k'-^e'T.^kZ-' (3.19) 

fc=0 

2 3 , , 

z = -C2 (3.20) 

r(3fc + ^) 

which agrees with the saddle point expansion. 
If we now take the logarithm we obtain 

log$(C) = ^log^-ilogC + ^d+ci^r^ + (c2-^c?)^C-' + ... (3.22) 



FN{g2) = ^-^/(to) + ^logy-^log$ (3.23) 



Comparison with ( |l.3| 



and (1.15) allows us to identify the coefficients ah, namely 





1 


- ao 


= +3 


-ai 


= arbit 


-a2 


= 3 ci 


-as 


= 9(C2 



1 2^ 



Moreover the singular term 



implies ( 1.2C| ). 



(3.24) 



^logC = ^logf(92-5*)^ (3.25) 



The partition function diverges at the critical value Qc = 92 = ~j- We shall show in the 
subsequent section that adding a term to f{t) (|3.3D, (f 



^53^^ 53 >0, small (3.26) 



is the most elegant way to come around this problem. 

4 Deformation of an catastrophe into curves of A2 catastro- 
phes 



Let 



fit) = \ogt - t - ^g2t - ^gst"" (4.1) 
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where 53 is a free parameter but we assume still that 



nto) + 



(4.2) 



then 33 defines a curve of A2 catastrophes. This curve has two real analytic branches intersecting 
at the A3 catastrophe (Figs 1 and 2), where 



/'"(to) 







In fact 

/'(to) = /"(to) 
can be solved in the form (see eqns (78), (79)) 

1 



to,± 



92 



1± (1 + 352)= 



and 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



-27(73,± = 2 + 992 + 2(1 + 352)2 

These are the curves drawn in Figs. 1 and 2 respectively. 

The partition function is convergent with (4J) if 53 > 0. We mentioned already in Section 3 
that an arbitrary small positive 53 can be used to give a well defined meaning to the t„-branch 
(at the dot). If we move from this point towards the j43 catastrophe we have 



< 2753 < 1 



< 52 < 



(4.7) 
(4.8) 



and remain continuously connected to the cuspoidal A2 catastrophe. This makes sense if we can 
steer the parameters 52 ; 53 at will. What if the system is such that it can adjust the parameter 
53 freely for fixed 52 ? Then it could jump (first order transition) to the t+-branch eventually. 
However, since at §2 = —j 



/(to,+) - /(to,-) 



92= 



log 3 



-0.2347) 



(4.9) 



the t_-branch is stable. 
We consider now 52 and 



as deformation parameters. Then 



'■o,± 



A 



-2/'(to 



A. 



+ 352 



(4.10) 



(4.11) 
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Our formulas are applicable in the case 



and 



52 < 0, A < 2 



(l--) +3(72 > 



Setting ( p^ ) equal zero gives the ^3 catastrophe since 



/"'(t0,±) = +2tn| 



1 



A\2 



+ 3^2 



(4.12) 
(4.13) 

(4.14) 



In order to avoid oscillating Airy functions which would lead to complex free energies we must 
moreover have 



+1 



Now we scale 



m 



sign{A/"'(to,±)} 

t — tQ = Xt] 

fit) = fito)-lA{t-to) + lf"'ito){t-tof + o{{t-tor 



so that 



or 



Nl 
~2Q' 



A 



1 

3 



Mf"'{to)\. 

The contribution of either branch to the partition function is (leading term only) 



2 N (92, 93, 



N_ 

IT 2 



± 



2 f -1 „i / 4 

P^^exp|y/(to,±)|to4 



N\f"'{to,±)\ 



where 



and 



C± = ^A^i|A| 



f"'ito,±) 



$±(C±) = dr/ exp{ ±C??+ ^??^} 



(4.15) 
(4.16) 
(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 
(4.22) 



C± are the contours obtained by a minimal deformation of the positively orientated real axis 
that makes the integrals convergent. We obtain 



$±(0 = vrBi(C) 
in either case. The constraint ( [4.15| ) renders C positive in either case. 



(4.23) 
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5 The general Am (cuspoidal) catastrophe 

We consider now the case 

gi = I, gk = 0, k>m + l (5.1) 

In this case we have 

= (-1)-^ ^ yt (5.2) 

The Am catastrophe occurs if 

= Vl<A:<m (5.3) 

This leads to the equations 



fk 

^0 l=k 



m 



This system can be solved in an elementary fashion for to and {gi} 
First we set 

(_l)ni-l 

k = m: — = Qm (5.5) 

which entails 

signg^ = (-1)™-! (5.6) 

since Iq must be positive 

io = l^ml"™ (5.7) 
Convergence of the partition function necessitates 

signc/m = +1 (5.8) 



which is compatible with (5^) only if m is odd. For even m we will employ the procedure 



discussed in the preceding section: We add a term 

-^gm+lt^^\ ff^+l\0 (5.9) 

to the action f{t). 

An intermediary step in solving the system (|5.4| ) is 

= (-I)'"-' Wc'^m, (l<A;<m) (5.10) 
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which for k = m is trivial. For A: = 1 we make use of = 1 to obtain to = In proving ( |5.1C 
we need the identity 



l-k 



l=k 



Denote the l.h.s. of (KU) by Then 



1 



which by Pascal's identity gives 



E(-i) 



l-k 



m 



/ m ' 

\kj 

^mk 



j_u fm — k^ 



Since P^h^^ = 1 ( 5.11 ) follows. Inserting ( ^.5[ ) into ( |5.10| ) we have finally 

5fc = (-1)'-' 
Now we deform this catastrophe by 



m 

, k . 



m 



\k-l 



m 

. k , 



9k = (-1, 
to = m + TQ 
Referring to translational invariance in t we postulate 

/(™)(to) = 

so that 



m ^ + Tk, {2 <k <m) 



fit) 



x:'^^^/«(to) + ^^"^"1' /(-+^)(to) + o((t-tor+^; 

^ /! (m + lj! V / 



We expand /'-'^(to) linearly in all to and r^, 2 < k < m. 
First we notice that 



0, < Z < m - 1 



ro=rfc=OVfc 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 
(5.16) 

(5.17) 
(5.18) 



(5.19) 



since 



(m) = 



(5.20) 
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k-r 



from (5.3). So there remains (ti = 0) 

f'^^\m + To) = -(^ - 1)! X! [ / 'm^'^Tk + quadratic terms in {r} 
Next comes the scaUng procedure 



k=l 



SO that 



leading to 



AT A^+i 



2 (m + 1)! 



f — io = ^f] 



(m) 



1 



m+1 



9 1 

/ ^ \ m+l 



Then introduce the scahng variables 



NX 



Ci = -jjf'^'Hm + To), (l</<m-l) 
which are kept 0(1) at the transition point by definition. It follows 



/(')(m + To) = Ai 



where 



and 



I 



m + l 



Ai = —^Ci = 0{l) 



so we have simple scaling of the derivatives implying that the coupling constants 

T2, T3, . . . , Tm-l, Tm 

scale along algebraic curves with ^ as single parameter, whereas tq is coupled to t„ 

To = (-l)'"m-r^ 
as follows from ( 5.17] ). In order to solve ( 5.21| ) for the Tk we denote 



Ik 



I - 1 



m 



k-l 
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Then the matrix inverse is 
It follows 



'l-V 
k-1 



m — 1 
,k-l , 



m 



1 < / < m - 1 
2< k<m 



(5.32) 



m—l 



Tk 



1=1 



_i Ai /iV\-fTi 



kl 



(/-I)! V2 



(5.33) 



Note that ti vanishes automatically. 

The contribution of this saddle point to the partition function is 



N_ 

vr 2 



r(f) 



^exp{^/(m)}(|^)'"+^$(Ci,C2,...,C 



2 



m—l I 



ZN{g) 

as leading term where 

/( Ti— 1 m+1 
dr/exp 5:a^?' + (-l)"^^ 



(5.34) 



'm—l 



(5.35) 



fc=i 



For odd m we identify C*-™^ with the positive real axis. For even m we define 



-Ci + - C 1 + C -1 

2 2 1 2(m+l) 2(m + l) 



(5.36) 



What are the conditions for asymptotic behaviour of ^ to be nonoscillating for large {Cfc} ? 
For m odd there is no problem since 



m odd : ^{C) > 



(5.37) 



by our definition. For even m only the case m = 2 has been dealt with already, e.g. by condition 
( 4.15| ) that leads to the correlation of signs in ( [4.22 ). For m > 4 studying the contributions of 
all subordinate catastrophes 

An , n < m 

is a complicated algebraic task. An approach giving an insight into the case m = 4 is presented 
in the Appendix. 

6 Differential equations and the renormalization equation 

The function <I>(Ci, C2, ■ ■ ■ , Cm-i) ( |5.35| ) satisfies the following set of linear differential equations 



Wk 



k £ {l,2,...m- 1} 



(6.1) 
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and 



(6.2) 



k=l 



In turn this system of m — 1 equations determines ^ to lie in the m-dimensional space of integrals 
( 5.34| ) with admissable contours C^"^\ 
Let us denote 

F(C) = logcI>(C) (6.3) 
Then -F(C) satisfies a system of nonlinear differential equations, e.g. for m = 2 



F" + {¥■ 



i\2 



c 



(6.4) 



which is the Airy differential equation in logarithmic camouflage. 

Independently of these differential equations F satisfies the renormalization group equation 







(6.5) 



k=2 



where we used the deviations of the coupling constants gk off their critical values (|5.15| ). 
Knowledge of the variables {Qi} as functions of and {rt} allows us to determine /^^(t) (insert 
(|5^ ) and (|]2|) into ( F^ ) 



/3fe(r) 



r=0 



We find 



N-^F = y aiC,i^F 



■dCi 



±F = y^±F 

dTk V ^"^fc ^0 



(6.6) 

(6.7) 
(6.8) 



and from 



, fH) 



ilk k ) OL,i 

Setting each coefficient of equal to zero gives 

(yiJ^lk = J^lrirk 
r 

with the susceptibility matrix 



5/3. (r) 



dTk 



Tl=0\/l 



(6.9) 



(6.10) 



(6.11) 
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So the matrix Nir diagonalizes this susceptibility matrix and {cr;} are its eigenvalues. In our 
inverse approach we can calculate by 

irk = Y.K^'^l^lk (6.12) 

I 

The sum can be performed and gives 

f3k{T) = -iM, J^TT2 + —-Tk + ^-{l-5km)Tk+l + 0{T^) 6.13 

\k — \ I m + 1 m + 1 m + 1 

{2<k<m) 



Appendix: Generalized Airy functions for Am, m even 

The asymptotic expansions of generalized Airy functions for large arguments as defined in ( ^.35| ) 
are treated with the same saddle point techniques which produced them. An Airy function for 
a catastrophe Am obtains contributions of all catastrophes j4/, \ < I < m. 
We restrict the arguments by 

Ci = 0, k <i<m-l (A.l) 

and consider corresponding reduced phase functions 

fkiv) = E0^' + ^^'"+' (A.2) 
^ m+1 

We let 

Ci = elCil, ICil^oo, = 1 (A.3) 

and couple the remaining variables to |Ci| 

Q = a\Ci\—^, 2<i<k (A.4) 
?7o = wiCil" (A.5) 

where tjq is the position of the Ai catastrophe, so that 

{a2,a3,...,uj} G 

is fixed in the limiting procedure. Then 

fkivo) = '^o{a2,a3,---,^)\Ci\^ (A.6) 
k ^ 

^-0(02, as,..., w) = €LO + y^aiJ + — (A.7) 

1^2 "^ + 1 
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We introduce the shorthand 

^i(a2,a3, • • • ,w) = ^^'0(02,03, • • • ,w) 
Then the Ai catastrophes appear at 

'J'i(a2,a3,---,w) = 

In the case k = 1 this is 

^-1(0;) = e + io"" = 
and has m sokitions uj. Only that solution is relevant in our context for which 



(Ai 



(Ai 



(A.IO) 



^-1=0 

is maximal. For a real nonoscillating asymptotic expansion we need moreover 

^'o real, ^'2 < 

For k = 1 this is possible only if 

e = —1, u; = — 1 



(A.ll) 



(A.12) 



Let us consider the case A; = 2 in more detail and concentrate on the case (|A.11| ). Then there 
are two possibilitues 



e = —1, —00 < oj < 0, a2 



2uj 



or 



1 1 + 

e = +1, < UJ < (m — 1)^™, 02 = 

2u> 



(A.13) 



(A.14) 



The case ( |A.12| ) is contained in (|A.13| ) but ( |A.14D is a new branch. At the end of this branch 
we have 



to = {m — ly 

and an A2 catastrophe. Denoting 

Xk = 

the Ai catastrophes contribute as leading term 

$(Ci,C2,...,Cfc,0,...) eXolCil"^ 



'l'i=0 



(A.15) 



(A.16) 



2tt 



\X2\ ICi 



(A.17) 
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Now to the A2 catastrophes. They show up first at fc = 2 as we just found, and are located at 

= 0, (A.18) 



For k = 2 these conditions are solved by 







(m - 1) u''' 
1 



02 



-muj 



,m—l 



(A.19) 

(A.20) 
(A.21) 



We deform the condition (A.18|) by a new parameter and use the standard scaling technique 



^1 = 
e = 

but (|A.19| ), (|A.21| ) is maintained. With 

A = 

and 



Ai 

(m- l)a;™ + Ai 



Ixsl iCi 



z = signxsAi A|Ci| 
we obtain the leading term of the asymptotic expansion {k > 2) 



$(Ci,C2,...,Cfc,o,...) 



m+1 



Ixsl ICi 



(A.22) 
(A.23) 



(A.24) 



(A.25) 



(A.26) 



if Ai approaches zero in such a fashion that z ( |A.25| ) remains fixed. Moreover Xk is now obtained 
from ^'fc by restriction to ( A.18 ) and ( |A.19| ). 

It is obvious that this procedure can be extended to all Ai and k. 
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